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In this report I summarize contributions on the highest energy cosmic rays 
to the 28th International Cosmic Ray Conference in Tsukuba, Japan, Sections 
HE 1.3, 1.4 and 1.5 involved over 80 oral presentations and about 120 posters. 
This large body of work attests to the vitality of a field committed to resolving 
the mystery of the origin of ultra-high energy cosmic rays. These cosmic rays with 
energies well above 10 18 eV are messengers of an unknown extremely high-energy 
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O ■ universe. 
Q-f 1. Introduction 

6 , 

+3 . Ultra-high energy cosmic rays are the highest energy messengers of the 

present universe. The highest energy cosmic photons observed thus far only reach 
about 10 13 eV (see, e.g., [1]). Photons of higher energies loose a significant fraction 
of their energies due to pair production in the cosmic background radiation as 
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5h ' they traverse large regions of intergalactic space. In contrast, cosmic rays are 

observed with energies as high as 3 x 10 20 eV and with fluxes well above upper 
limits on high-energy photon fluxes (see Fig. 1). 

The origin of these relativistic particles remains a mystery hidden by the 
fact that cosmic rays do not point back to their sources. These charged parti- 
cles are deflected by magnetic fields that permeate interstellar and intergalactic 
space. Galactic magnetic fields are known to be around a few micro Gauss in 
the Galactic disk and are expected to decay exponentially away from the disk [2] . 
Intergalactic fields are observed in dense clusters of galaxies, but it is not clear 
if there are intergalactic magnetic fields in the Local Group or the Local Super- 
galactic Plane [3]. On larger scales, magnetic fields are known to be weaker than 
~ 10 nano Gauss [4]. As cosmic ray energies reach 10 20 eV per charged nucleon, 
Galactic and intergalactic magnetic fields cannot bend particle orbits significantly 
and pointing to cosmic ray sources becomes feasible (see, e.g., [5,6]). At ultra high 
energies there is finally an opportunity to begin cosmic ray astronomy. 

In addition to the ability to point back to the source position, cosmic ray 
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Fig. 1. Cosmic Backgrounds: photons in black [8] and cosmic rays in red circles. 



protons of energies around 10 20 eV should display a well-known spectral feature 
called the GZK cutoff [7]. This cutoff was proposed in 1966 by Greisen, Zatsepin 
and Kuzmin as a natural end to the cosmic ray spectrum due to photopion produc- 
tion off the then recently discovered cosmic microwave background radiation. The 
presence of microwave photons through cosmic space induces the formation and 
subsequent decay of the A + resonance for protons with energies above ~ 10 20 eV 
that traverse distances longer than ~ 50 Mpc. The effect of photopion production 
is to decrease the energy of protons from distant sources resulting in a harden- 
ing of the spectrum between 10 19 eV and 10 20 eV followed by a sharp softening 
past 10 20 eV. Depending on the maximum energy of ultra high-energy cosmic ray 
sources and their distribution in the universe, the spectrum may harden again 
past the GZK feature displaying the injected spectrum of nearby sources. 

The search for the origin of the highest energy particles is being undertaken 
by a number of experiments. At the 28th International Cosmic Ray Conference 
in Tsukuba, Japan, the two largest exposure experiments, the Akeno Giant Air- 
shower Array (AGASA) and the High Resolution Flys Eye (HiRes) presented a 
number of interesting results with additional talks by the Yakutsk array and a 
number of upcoming large experiments such as Auger and EUSO. AGASA is a 
100 km2 ground array of scintillator and muon detectors. Their collaboration 
presented 3 talks and a few posters. HiRes is composed of fluorescence telescopes 
built in two different sites in the Utah desert to be used as a stereo fluorescence 
detector. The longer exposure site, HiResl, has 21 mirrors while the second site, 
HiRes 2, has 42 mirrors. The HiRes collaboration presented 6 oral and 5 posters. 
Their main results are discussed in §2 as the present state of UHECR results 
are summarized. In §3, some implications and caution with the present data are 
presented. In §4 observatories now under construction are discussed. The Pierre 
Auger Project gave a preview of the powerful results to expect in the near future 
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as construction is completed. Their 12 oral and 17 poster presentations showed 
the potential power of a large hybrid detector of fluorescence and water tanks. 
Another large project now in the development phase includes the Telescope Array 
project. The Extreme Universe Space Observatory (EUSO) had a large presence 
in the 28th ICRC presenting 11 talks and 17 posters. With a fluorescence tele- 
scope in the International Space Station, EUSO will cover a very large volume 
during its operations (a base of 160,000 km 2 ) focusing on the highest energy re- 
gion (well above 5 x 10 19 eV). In §5 we conclude by looking forward to the 29th 
ICRC in 2005. 

2. Present State of UHECRs 

Observations of cosmic rays at the highest energies have yielded measure- 
ments of the spectrum, arrival direction distribution, and composition of UHE- 
CRs. We start by highlighting new results on the composition, then follow with a 
discussion of arrival directions, and conclude this section with the measurements 
of the spectrum. 

2.1. Composition 

A comprehensive study of the composition of UHECRs with energies be- 
tween ~ 10 18 eV and 10 19 eV was presented by the HiRes Collaboration using 
Stereo data [9]. They find that the elongation rate agrees best with a light com- 
ponent above 10 18 eV. The HiRes Stereo study fits well with the HiRes Prototype- 
MIA analyses that shows a trend from heavier to lighter nuclei as the energies 
grows from 10 17 eV to 10 18 eV. The elongation rate versus energy is reproduced 
here in Fig. 2. The composition studies do not reach energies near the GZK 
feature but should be available by the 29th ICRC. this trend indicates that a new 
light component, most likely extragalactic protons, become the dominant compo- 
nent at the highest energies in this range. The change of slope seen in the spectra 
(as discussed below) may also indicate the transition to extragalactic. 

The AGASA collaboration presented constraints on the UHECR compo- 
sition [10] from muon density studies. They also find that a lighter component 
may be dominant at higher energies. For energies between 10 19 eV and 10 19 3 eV 
the iron fraction is less than 35% while for energies above 10 19,3 eV, they expect 
less than 76% as iron. These findings are consistent with studies in the Akeno 
array [11], Haverah Park [12], Volcano Ranch [13], and HiRes [9]. AGASA also 
set upper limits on a gamma-ray component of less than 34% for energies above 
10 19 eV, and less than 56% for energies above 10 19 ' 3 eV [10]. 
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Fig. 2. Elongation rate versus energy from Stereo HiRes [9]. 

2.2. Distribution of Arrival Directions 

The arrival directions of cosmic rays are strongly isotropized by the effect 
of magnetic fields on the propagation of charged primaries. As primary energies 
increase, the isotropizing effect of cosmic fields should ease and the location of 
sources should be revealed. Thus, anisotropies on the distribution of arrival direc- 
tion hold the key to identifying sources of UHECRs. We should expect large-scale 
anisotropies associated with the large-scale source distribution and small scale 
clustering associated with point sources. 

The first hints of large-scale anisotropies have been reported for energies 
around 10 18 eV by AGASA [14] where a 4% excess over the diffuse flux effect is 
reported at the 4 a level. HiRes [15] and SUGAR [16] confirm an excess with 
lower significance. The anisotropy points towards the Galactic Center region, 
which indicates a Galactic UHECR component. This small anisotropy may be 
associated with cosmic ray neutrons that originated either in the Galactic Center 
region or in the Cygnus spiral arm. 

If protons are the dominant component of cosmic rays at the highest en- 
ergies (as indicated in $2.1), the arrival distribution should show anisotropies 
associated with the source distribution as the primary energy increases. Hints of 
small scale clustering at the highest energies were reported by AGASA [17]. The 
two point angular correlation for events above 4 x 10 19 eV and zenith angle less 
than 50 degrees shows a peak at small angles, which corresponds to 6 doublets 
and 1 triplet. (An analysis with looser cuts yields 2 triplets and 6 doublets with 
energies between 10 19 eV and 5 x 10 19 eV.) The ratio of events in clusters to all 
events seems to increase with energy as oc E' 1 '^ ' 5 . For energies above 5 x 10 19 
eV, where the GZK feature is expected, the statistics is too low for observing 
clusters. There is no significant arrival time correlation to indicate a source time 
dependence. 
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Fig. 3. AGASA and HiRes exposures (courtesy of HiRes and AGASA collaborations). 



In addition to studying the angular correlation of the events, AGASA has 
displayed the clustering in a two-dimensional distribution of right ascension and 
declination. They find directionality in this two-dimensional distribution that 
can be associated with distortions caused by the Galactic magnetic field. This 
polarization effect on the arrival distribution of point sources should be evident 
if clusters with larger statistics are observed by future observatories. The study 
of the distortion shape of point source will place interesting constraints on the 
structure of the Galactic and Intergalactic magnetic field [6] . 

In contrast to the clustering reported by AGASA, HiRes sees no significant 
departure from an isotropic distribution. The exposure of Stereo HiRes is not yet 
comparable to the AGASA exposure, but Monocular HiResl data has comparable 
exposures at the highest energies (see Fig. 3). The Monocular HiRes analysis is 
complicated by the elongated shapes of each event. A sophisticated analysis of 
small scale clustering of HiResl data controlled by simulations that takes the 
asymmetric sky exposure into account yielded no significant clustering with an 
upper limit of 4 doublets (at 90 % c.l.) [18]. They also report upper limits on 
point sources in Cygnus X-3 and on a dipole anisotropy in the direction of the 
Galactic Center, Centaurus A, and M-87 [18]. (Given the partial coverage of 
the sky, only a limited search for a dipole anisotropy can be conducted.) No 
significant clustering was seen also in the HiRes Stereo data, consistent with the 
Monocular finding [19]. 

If the AGASA clusters of events are the first hints of point sources, the 
number density of sources can be estimated to be around 10 -5 Mpc~ 3 [20] if point 
sources have injection spectrum oc E~ 2M . The same simulated point sources for 
the Auger exposure after 5 years would yield 70 clustered events with energies 
above 10 20 eV while EUSO would observe 180 to 360 clustered events in 3 years 
above 10 20 eV. A definitive measurement of point sources is well within the reach 
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Fig. 4. AGASA spectrum 

of the next generation observatories 
2. 3. Spectrum 

The cosmic ray spectrum past 10 20 eV will indicate the presence or absence 
of the GZK feature, which can be related to the type of primary (e.g., protons) 
and source (injection spectrum and spatial distribution) of UHECRs. The largest 
two experiments have conflicting findings past 10 20 eV. AGASA does not see a 
steepening of the flux above 10 20 eV as expected from a GZK cutoff and reports 
a continuation of the lower energy spectrum and 11 Super-GZK events, i.e., 11 
events with energies above 10 20 eV [21] (see Fig. 4). These findings argue against 
the notion of extragalalactic proton sources and for a new and unexpected source 
of UHECRs. 

In contrast, the HiResl Monocular spectrum indicates smaller fluxes past 
10 20 eV which is consistent with a GZK feature [22]. HiRes reports only two 
events with energies above 10 20 eV (1 seen in Stereo). Spectral index fits show 
that the HiRes spectrum is not fit well by a single power law (see Fig. 5). The 
systematic discrepancy between HiRes and AGASA spectra is ~ 30% in energy 
scales. Still unclear is also the energy calibration between HiRes and Flys Eye: 
HiRes energies are ~ 7% lower when compared to Flys Eye. 

HiRes collaboration also analyzed part of their Stereo data and showed a 
energy spectrum [23]. The significance of any feature in these data is still weak 
given the low number of events analyzed, however, a change in spectral index can 
be seen around 10 18 ' 6 eV (see Fig. 6), where a simple two component composition 
is folded in (76 % of proton and 24 % iron mixture). This change in spectral 
index may be due to the Galactic to extragalactic component transition. 

The implications of the differing results from AGASA and HiRes are espe- 
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Fig. 5. HiRes 1 and 2 Monocular spectra. 
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Fig. 6. HiRes Stereo spectrum. 



daily intriguing at the highest energies. AGASA reports no indication of a GZK 
feature with the observation of 11 Super-GZK events. In contrast, HiRes 1, which 
is the HiRes site with comparable exposures to AGASA at the highest energies, 
does not report observations of as many Super-GZK events. HiRes was designed 
to be used as a stereo detector in order to best control the determination of each 
shower axis, but the exposure of HiRes 2 and, thus, of Stereo HiRes is at present 
significantly below that of AGASA (see Fig. 3). 

Possible sources of systematic errors in the energy measurement of the 
AGASA experiment were comprehensively studied to be at around 18 % [24]. 
Systematic errors in HiRes are still being evaluated, but are likely to be dominated 
by uncertainties in the absolute fluorescence yield, the atmospheric corrections, 
and the calibration of the full detector. Laser probes were used to calibrate some 
of the atmospheric uncertainties at the HiRes sites [25]. 
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Fig. 7. AGASA statistics [26]. 
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Fig. 8. HiRes statistics [26]. 



Although control of systematic errors is crucial, the statistics accumulated 
by both HiRes and AGASA are not large enough for a clear measurement of the 
GZK feature. Figs. 7 and 8 show the range of 400 simulated spectra of protons 
propagating in intergalactic space with injection spectral index of 2.8 for AGASA 
(Fig. 7) and 2.6 for HiRes (Fig. 8). Uncertainties due to the stochastic nature of 
the photo-pion production dominates the spectrum around the GZK feature for 
experiments low statistics (see Figs. 7 and 8) [26]. The GZK feature can only be 
clearly measured with larger exposures above 10 20 eV. 

If systematic errors at HiRes are also around 20 % the two data sets are not 
totally inconsistent. Once the low statistics of both detectors together with 15 % 
systematics are taken into account, the disagreement between the two experiments 
is around 2 a [26] (see Fig. 9). The low exposure of both experiments prevents 
an accurate determination of the GZK feature or lack of it. 
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Fig. 9. AGASA with -15% energy shift and HiRes with +15% shift [26]. 
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Fig. 10. Yakutsk spectrum [27]. 



Finally, Yakutsk also presented a spectrum that is also somewhat dis- 
crepant with AGASA and HiRes, as seen in Fig. 10 [27]. 

3. Lessons for the Future 



The present state of UHECR data can be summarized as follows. Between 
10 18 eV and 10 19 ' 3 eV the composition shifts to lighter component and the dis- 
tribution of arrival directions is mainly isotropic. AGASA sees evidence of small 
scale clustering that in not seen at HiRes. The spectra from AGASA and HiRes 
show a systematic shift through the range of observed energies and the statistics 
at the highest energies is not sufficient to determine the existence of a GZK cutoff. 
When the two spectra are plotted on a flux versus energy plot (see Fig. 11), the 
discrepancies are not as accentuated as in the traditional plots of flux times E 3 . 
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Fig. 11. AG AS A and HiRes spectra (thanks to D. Bergman). 



The lessons for the future are clear: improve the statistics significantly at the 
highest energies and understand the sources of systematic errors. 

3.1. Fluorescence in the Lab 

The results from AGASA and HiRes show the need to understand and 
control systematic effects within each technique and to cross-calibrate the two 
techniques presently available for ultra-high energy cosmic ray (UHECR) studies. 
One step towards the control of systematic effects is the study of fluorescence 
yields in the laboratory. Following on previous measurements of the total yield of 
fluorescence lines (between 300 and 400 nm [28]), N. Sakaki et al. [29] measured 
individual spectral lines in six wavelengths and estimated a systematic error of ~ 
13.2% in overall yield and in individual lines. 

In addition to measuring fluorescence yields, laboratory experiments are 
now being developed to study airshowers in the laboratory. The Fluorescence from 
Air in Showers project or FLASH is using the SLAC beam to produce showers in 
a controlled experimental setting [30]. They presented preliminary results from 
a thin target run and expect to reach a lower than 10% accuracy in the total 
fluorescence yield and individual spectral lines. 

3.2. Model Uncertainties 

One of the main sources of uncertainty in ground array energy reconstruc- 
tion is fits to airshower simulations. A number of presentations showed the im- 
provements in speed and accurate physics input of the new generation of airshower 
simulators. The importance of using an accurate description at low energies, such 
as provided by FLUKA, for water tanks far from the shower core was highlighted 
[31], while hybrid codes were shown to improve performance speed [32]. However, 
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Fig. 12. Elongation rate on a range of QGSJet models [33]. 
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Fig. 13. Showers in atmospheres in different seasons [34]. 

the strongest message was that different experiments should use the same simu- 
lation codes in their data analysis such that different shower simulators do not 
introduce added systematic effects. 

The dependence on shower simulations is not only relevant to energy cali- 
bration in ground arrays but it is also crucial for distinguishing between light and 
heavy baryonic compositions of the primaries. An extreme comparison in com- 
position studies using QGSJet shows the range in model prediction of elongation 
rates as in Fig. 12 [33]. 

In addition to biases due to simulations of particle interactions, composi- 
tion studies may also suffer from porr modeling of the atmosphere. In particular, 
the use of "Standard" atmospheres can bias composition studies [34]. As see in 
Fig. 13, iron in the winter may look like protons in the summer. 

Uncertainties in interpreting cosmic ray data may also arise from astro- 
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Fig. 14. GZK feature B dependence [35]. 



physical unknowns. The shape of the GZK cut-off depends the source distribution 
and the magnetic fields in the Galaxy and in the extragalactic medium. Fig. 14 
shows variations on the predicted GZK feature due to different magnetic field 
structures [35]. 

The presence of a feature depends on the sources being astrophysical ac- 
celerators versus new physics models of decay from high energies. If sources are 
astrophysical and the primaries are protons, the photo-pion production, which 
generates the GZK feature, will signal a new source of ultra-high energy neutri- 
nos. These cosmogenic neutrinos can be used to study fundamental interactions 
at center of mass energies well above a TeV. In particular, the existence of large 
extra-dimensions above a TeV scale will modify neutrino generated showers [36] . 

4. Preview of the Next Generation 

Neither AGASA nor HiRes have the necessary statistics and control of 
systematics to determine in a definitive way the existence of either the GZK 
feature or of a novel source of Super-GZK events. Moreover, if the AGASA 
clusters are an indication of point sources of UHECRs, a large number of events 
per source will be necessary to study their nature. In order to discover the origin 
of UHECRs, large new projects are now under construction including the Pierre 
Auger Project, the Extreme Universe Space Observatory, and the Telescope Array. 
In addition, projects such as CHICOS, SCROD, and ASHRA also presented lower 
budget alternatives that may contribute to the field. 

4-1. Pierre Auger Observatory 

The Pierre Auger Project will consist of two giant airshower arrays one 
in the South and one in the North each with 1600 water Cherenkov detectors 
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Fig. 15. Auger South statistics at the GZK feature [26]. 



covering 3000 km 2 and four sites of fluorescence telescopes. Auger is being built 
to determine the spectrum, arrival direction, and composition of UHECR in a 
full sky survey. The survey should provide large event statistics and control of 
systematics through detailed detector calibration of the surface array and flu- 
orescence detectors individually in addition to the cross-calibration of the two 
detection techniques through the observation of hybrid and stereo-hybrid events. 
Depending on the UHECR spectrum, Auger should measure the energy, direction 
and composition of about 60 events per year above 10 20 eV and about 6000 events 
per year above 10 19 eV. In addition, it should be able to detect a few neutrino 
events per year if the UHECRs are extragalactic protons. 

The Auger collaboration consists of about 250 scientists from 16 countries. 
The collaboration presented 12 oral papers, including a highlight talk, and 17 
posters that described the progress in the construction of the Southern observatory 
in the province of Mendoza in Argentina and some findings from the prototype 
Engineering Array with about 40 tanks and two fluorescence telescopes. At the 
time of this meeting there were about 140 tanks deployed in the array [36b]. 

The Auger surface array is composed of stand alone 1.5 meter tall water 
tanks that are powered by solar cells, timed by GPS systems, and communicate 
via radio antennas. Three photomultipliers per tank register the Cherenkov light 
when shower particles cross the tanks. Having three photomultipliers per tank 
allows the self-calibration of each tank in the field [37]. The height of the tanks 
makes the ground array an excellent detector for inclined showers. Inclined show- 
ers and their asymmetries [38] allow for a novel method for composition studies 



Analysis of the Engineering Array surface detector events fits well with 
the shower lateral distribution from Monte Carlo simulations [40] and with the 
fluorescence detector analysis of hybrid events [41]. Hybrid detection is a powerful 
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measurement of individual showers and can be used to reach large statistics on 
energies down to 10 18 eV with the use of fluorescence and a small number of tanks 
per event [42]. The ability to study events at 10 18 eV in the Southern hemisphere 
will be crucial in confirming the reported anisotropics toward the Galactic Center 
region. The combination of mono fluorescence events that triggered even a single 
tank allows for great angular reconstruction of events comparable to stereo events 
[43]. 

The fluorescence detectors at the Auger observatory have a complete cali- 
bration system. The atmospheric monitoring includes lasers, lidars, ballon radio 
sondes, cloud monitors, and movable calibration light sources [44,45]. In addi- 
tion, the whole telescopes including mirrors are calibrated from front to end with 
light sources. The analysis procedure for the fluorescence telescopes include an 
iterative process that can model the direct Cherenkov contribution [36b]. With 
all these technical developments Auger should be able to show some interesting 
results by the next ICRC meeting. 

4-2. Telescope Array 

Another upcoming experiment is the recently approved Telescope Array 
(TA) which presented 3 oral contributions and 1 poster. TA's plan is also to 
have a hybrid detector of three fluorescence telescopes overlooking a scintillator 
array. The array would cover about 400 km 2 with an array with 1.2 km spacing 
to concentrate their studies in the transition between Galactic to extragalactic 
UHECRs from below ~ 10 17 eV to just above ~ 10 20 eV. TA should be able 
to see some super-GZK events but with significantly smaller statistics than the 
Auger project. Instead, TA is planning to concentrate their efforts in having a 
broad reach in energies to study the spectrum during the transition from Galactic 
to extragalactic and heavy to light primaries. 

4.3. EUSO 

The Extreme Universe Space Observatory (EUSO) is a fluorescence de- 
tector designed for the International Space Station (ISS) aiming at observations 
of extremely high energy cosmic-rays (EHECRs), i.e., cosmic rays between 10 20 
and 10 22 eV. EUSO will see showers from abovethe atmosphere and will have full 
sky coverage due to the ISS orbit. This project is a good complement to ground 
arrays since it will focus on larger energy scales and will have different systematic 
effects. Their threshold may be above 5 x 10 19 eV depending on technical features 
of the fluorescence detectors. 

The EUSO collaboration includes about 120 scientists from 7 countries. 
The collaboration presented 10 oral papers and 18 posters. They are now in 
phase B at NASA and ESA. The telescope's expected angular resolution is about 
0.2 degrees and the energy resolution about 14%. The aperture may reach 3 x 10 6 
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Fig. 16. EUSO statistics at the GZK feature [26] . 



km 2 -sterad with a 10% duty cycle. This can translate into about 3000 events per 
year for energies above 10 20 eV (see Fig. 16). For a source density of 10 -5 Mpc~ 3 
EUSO would observe one event per source in half a year versus about 6 years 
with Auger South [20]. The 200 Giga-ton atmosphere is also a good neutrino 
detector. Horizontal, Earth skimming and upward going tau neutrinos are some 
of the possible ways neutrinos could be seen by EUSO. The challenge for EUSO 
with neutrinos is the low energy threshold since cosmogenic neutrinos have typical 
energies well below 10 19 eV. The possibility of detecting direct Cherenkov light 
from upward going showers opens a new possibility in neutrino detection with 
EUSO. 

5. Conclusion 

After decades of attempts to discover the origin of ultra-high energy cosmic 
rays, present results are still inconclusive. The two largest exposure experiments, 
the Akeno Giant Airshower Array (AGASA) [10] and the High Resolution Flys 
Eye (HiRes) [11], show a systematic off-set in the observed spectrum indicating 
systematic errors in either the energy or the flux determination of one or both 
experiments (see Fig. 11). 

The results from these experiments show the need to understand and con- 
trol systematic effects within each technique and to cross-calibrate the two tech- 
niques presently available for ultra-high energy cosmic ray (UHECR) studies. In 
addition, the lack of sufficient statistics limits the discussion of an excess flux or 
a drop in flux around the GZK feature. Next generation experiments are gearing 
up to accumulate the necessary statistics while having a better handle on the 
systematics. In the following decade, we may see the growth of a new astronomy 
with ultra-high energy charged particles. 
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